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Abstract. Inspired by Axelrod's model of culture dissemination, we introduce and analyze a model for a 
population of coupled oscillators where different levels of synchronization can be assimilated to different 
degrees of cultural organization. The state of each oscillator is represented by a set of phases, and the 
interaction -which occurs between homologous phases- is weighted by a decreasing function of the distance 
between individual states. Both ordered arrays and random networks are considered. We find that the 
transition between synchronization and incoherent behaviour is mediated by a clustering regime with rich 
organizational structure, where some of the phases of a given oscillator can be synchronized to a certain 
cluster, while its other phases are synchronized to different clusters. 

PACS. 05.45.Xt Synchronization; coupled oscillators - 87.23.Ge Dynamics of social systems - 89.75.Fb 
Structures and organization in complex systems 



1 Introduction 

Robert Axelrod's model for the dissemination of culture 
PQ aims at giving a simplified picture of the processes that 
shape the distribution of cultural features in a society, in 
particular, the balance between the trend to convergence 
due to cultural affinity, and the mechanisms that main- 
tain diversity. This model has attracted the attention of 
physicists because of its nontrivial phenomenology, which 
includes the existence of a variety of absorbing states, crit- 
ical phenomena, and non-monotonic dynamics [2j. Among 
other topics, the transition between homogeneous and het- 
erogeneous from states under variation of the number of 
cultural traits [3], the effects of noise [4] and non-trivial 
underlying topologies [5], as well as the mean-field limit 
[H] have been characterized. 

In its original formulation pQ, Axelrod's model de- 
scribed the cultural profile of each agent in a population 
by a vector of features. Each feature can adopt a num- 
ber of traits. Agents are situated at the nodes of a bidi- 
mensional square array, and interaction occurs between 
nearest neighbours. Two neighbour agents interact with 
a probability proportional to their cultural overlap, i.e. 
to number of features in their profiles whose traits coin- 
cide. As the result of the interaction, the profiles change, 
approaching each other. If the two profiles exhibit com- 
pletely different traits, interaction is impossible. This sys- 
tem reaches a frozen state where any two neighbours have 
either identical or completely different profiles. The final 
state can be homogeneous, with all the population having 
the same profile, or consist of coexisting domains, thus 



maintaining cultural diversity [1,2J. It turns out that, as 
the number of possible traits per feature grows, there is a 
transition from homogeneity to diversity [3]. 

From a physicist's viewpoint, there is a strong similar- 
ity between the self-organization process that determines 
the distribution of cultural profiles in Axelrod's model, 
and the evolution of a spatially extended system able to 
develop either homogeneous states or phase coexistence. 
Synchronization of distributed interacting oscillators is an 
example (TBS]. In these systems, the attractive interac- 
tions that may lead to the formation of synchronized do- 
mains plays a role similar to the tendency of interacting 
individuals to become culturally more similar in Axelrod's 
model. On the other hand, diversity between individual at- 
tributes of oscillators impedes their synchronization, like 
too much cultural divergence does not allow interaction. 

Inspired by Axelrod's model, in this paper we intro- 
duce and analyze a system of interacting oscillators, where 
the state of each oscillator is characterized by a set of sev- 
eral phases. Each phase would correspond, in Axelrod's 
model, to a cultural feature. Homologous phases of differ- 
ent oscillators are coupled to each other, with a coupling 
intensity that depends on the overall distance between the 
states of the two oscillators. We show that, in this system, 
the synchronization transition can be straightforwardly 
identified with Axelrod's transition from cultural homo- 
geneity to diversity. Additionally, we disclose an inter- 
mediate regime with rich organizational structure, where 
clusters of mutually synchronized oscillators, different for 
each phase, spontaneously appear. 
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2 Model 

We consider a population of N oscillators, each of them 
occupying a node in a network. The state of oscillator i is 
characterized by a set of F phases, 4>{(t) £ (0, 2w) (/ = 
1, . . . ,F), The dynamics of the phases 4>{(t) is governed 
by Kuramoto-like equations [TUp] 
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(4-4l 
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where the sum runs over the oscillators j connected to i 
by network links, which define the neighbourhood A/j of 
i. The number of oscillators in A/i is Vi. The non-negative 
coupling constant kij is computed as a prescribed function 
of the distance D^ between i and j, 



k l3 = IMI),,). 
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The distance Dij , in turn, characterizes the difference be- 
tween the individual states of oscillators i and j, as 



Di 



1 



F 

2F E 



1 — cos( 



(3) 



Note that D^ is non-negative, and = if and only 



if 



, F. Moreover, the maximal 



j = 4>{ for all / = 1, 
possible value of the distance is £>y = 1 . For the function 
K(D), which defines the coupling constant through Eq. 
([2]), we choose 

1 - (aD) r for D <a~\ 
K(D) = { (4) 
otherwise, 

with a > 1 and r > 0. Thus, the coupling constant is 
maximal, kij = 1, when the distance Dy vanishes, and 
decreases monotonically as Dij grows. It reaches zero for 
= cv _1 , and kij — for larger distances. The exponent 
r controls the behaviour of kij at small distances: for r < 1 
and r > 1, the coupling constant at Dy = displays, 
respectively, a cusp and a flat maximum. For r = 1, the 
dependence of kij on the distance is linear. 

When the coupling constant kij is positive, the sum- 
mand in the right-hand side of Eq. ^ represents an at- 
tractive interaction between 4>{ and 4>j. Under its action, 
the two phases will tend to become mutually synchronized. 
While the functional dependence of the interaction couples 
phases separately for each value of / = 1, . . . , F, its inten- 
sity -determined by the coupling constant- is given by the 
joint contribution of the difference of all the F phases be- 
tween the two oscillators through their distance Dij, Eq. 

The multi-phase oscillator system defined by Eqs. (U 
to (|4]) can be qualitatively compared to Axelrod's model as 
follows. The F phases of each oscillator are associated with 
an agent's cultural features. Different values of each phase 
-which, in our system, vary continuously over (0, 2n)- cor- 
respond to the traits of Axelrod's model. The distance Dij 



gives a measure of the cultural dissimilarity between two 
agents. If it is too large, the two agents do not interact. 
For small distances, on the other hand, the intensity - 
or, in Axelrod's model, the probability- of the interaction 
grows as the distance decreases. Synchronization between 
phases of two or more oscillators can be assimilated to 
their cultural consensus. This may occur with respect to 
all the phases -i.e., all the cultural features- or just some 
of them. A given oscillator could have some of the phases 
synchronized with a certain part of the population, and 
the other phases synchronized with another part, thus giv- 
ing rise to a state of partial synchronization equivalent to 
a rich diversity of cultural domains. As in other variants of 
the model [2J, we are here considering that the structure of 
the population is defined by a generic network, instead of 
the bidimensional array originally considered by Axelrod. 

In our system, we expect that collective synchroniza- 
tion is easier to achieve when the interaction range with 
respect to the distance Dij between the oscillators' states, 
given by a -1 in Eq. Q, is larger. This corresponds, in 
Axelrod's picture, to a small number of traits per cultural 
feature. It has been shown, in fact, that Axelrod's model 
exhibits a transition between cultural homogeneity and 
diversity as the number of traits in each feature grows [HI 
[2] . Similarly, in Eq. |(TJ) a synchronization transition is ex- 
pected to occur depending on the interaction range oT 1 
being larger or smaller than the typical distance between 
the states of any two oscillators in the population. 

Assume that the system is in a fully unsynchronized 
state, with all the F phases of all oscillators uniformly 
distributed over (0,27r). The expected average value for 
the distance Dij between any two oscillators is 
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If this mean distance is larger than a -1 , the unsynchro- 
nized state should be stable as oscillator pairs do not inter- 
act on the average. Consequently, we predict a transition 
from synchronization to desynchronization at the critical 
point 

a c = 2. (6) 

The width of this transition should be controlled by the 
dispersion in the values of D^ . Under the same assump- 
tions as for the calculation of (Dij), the mean square dis- 
persion of the distance turns out to be 



°n = (Dij) 2 ] 
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According to this estimation, thus, the transition width 
should decrease as the number of phases F grows. On the 
other hand, since our estimation is based on the evaluation 
of the average distance between just two oscillators, the 
width does not depend on the population size N. 

In the following section, we present numerical results 
for our model, Eqs. {J) to (@]), for populations distributed 
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over ordered and random networks. We confirm the above 
predictions on the synchronization transition, and show 
that around the critical point the population segregates 
into clusters of mutually synchronized oscillators, display- 
ing a high degree of organizational diversity. 



3 Numerical results 

The results presented in this section are the outcome of nu- 
merical calculations made on populations of various sizes, 
ranging from N — 100 to 1000 oscillators. As for the un- 
derlying interaction network, we consider a one-dimen- 
sional ordered array and random Erdos-Renyi networks 
1 1 1 j . In the former, each node is connected to its four near- 
est neighbours, two to each side. In the latter, the average 
number of neighbours per node equals four. 

Bearing in mind the analogy with Axelrod's model, we 
focus our attention on the organization of the population 
into synchronized clusters at asymptotically long times. A 
synchronized cluster is a group of oscillators whose phases 
coincide. Given that we deal with multi-phase oscillators, 
where the state of each of them is described by F phases, 
we analyze the organization into groups for each phase 
separately. Numerically, we consider that two phases are 
synchronized if they differ in less than A<p = 10~ 4 . To 
statistically characterize the organization into groups we 
measure the normalized mean number of groups g and the 
dispersion a g , averaging over phases and realizations. The 
normalization of the number of groups is performed with 
respect to its maximum value, N. 

A more detailed characterization of the organization 
into groups is achieved by introduced a normalized Ham- 
ming distance, h, averaged over phases and oscillator pairs, 
as follows. First, we represent the synchronization state of 
each oscillator pair by defining a set of F matrices M* 
(/ = 1, . . . ,F), N x N in size, with elements 

!1 if <j){ and M are synchronized, 
(8) 
otherwise. 

The Hamming distance between two of these matrices, 
and , is defined as 

N 

H f f>= £ |m£-m£|. (9) 

It equals zero if and only if the two matrices are identical, 
and its maximum value is N(N—1), The normalized Ham- 
ming distance h is obtained by averaging Hfp over all the 
matrix pairs and normalizing to the maximum value. 

Figure Q] shows results for the normalized number of 
groups g as a function of the parameter a of Eq. |(4]) , for 
various values of the exponent r and both network topolo- 
gies. A feature common to all curves is the rather sharp 
transition from small to large values of g as a grows, just 
above the critical value predicted in the preceding section, 




a 



Fig. 1. Normalized number of groups as a function of a for the 
one-dimensional array (grey) and for random networks (black) 
with N = 100 and F = 50. Each set of curves corresponds to 
three values of the exponent r = 0.1 (solid), 1 (dashed), and 3 
(dotted). 



Eq. (J6j) - This corresponds, as expected, to a desynchro- 
nization transition as the range of coupling decreases. To 
the right of the transition, in fact, the number of groups 
equals the population size, reveling a state where all os- 
cillators are mutually unsynchronized. 

Around and above the transition the behaviour is rather 
independent of the parameters and of the underlying topol- 
ogy, but a substantial difference is apparent for a < 2. 
While in random networks the normalized number of groups 
remains small and close to its minimal value g = 1/N, 
revealing full synchronization of all oscillators in all their 
phases, for the one-dimensional array g attains rather large 
values, g > 0.5. By inspecting single realizations on the 
one-dimensional array, we have verified that these large 
values of g are due to the occasional occurrence of the so- 
called twisted states [12]. In twisted states, at asymptoti- 
cally long times, oscillator phases are not synchronized but 
vary linearly along the array. These spatially-correlated 
distributions of phases are possible due to the underlying 
topological order. From the viewpoint of our approach, 
however, they represent fully unsynchronized states where 
the number of groups equals the size of the population, N. 
Therefore, their contribution makes the average normal- 
ized number of groups g grow. 

For larger values of a, as illustrated by FigfH differ- 
ences between results for random networks and the one- 
dimensional array, as well as for different values of r, are 
much less pronounced. Thus, in our analysis for a > 2, we 
focus attention on random networks and fix r = 1. 

We first verify the prediction made in the preceding 
section, Eq. J7]), that the width of the transition should 
decrease with as the number of phases F grows. This is 
confirmed by the numerical results shown in Fig. [21 where 
we plot the normalized number of groups g as a function of 
a for three values of F. Figure [3] shows that, in agreement 
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Fig. 2. Normalized number of groups g as a function of a for 
TV = 100 and F — 50 (solid), 10 (dashed), and 5 (dotted). 



Fig. 4. Histograms of the normalized number of groups g for 
several values of a, corresponding to the realizations for TV = 
100 used to calculate the averages and mean square dispersions 
shown in Fig. [3] 
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Fig. 3. Normalized number of groups g as a function of a for 
F = 50 and TV = 100 (solid), TV = 300 (dashed), and TV = 1000 
(dotted). Empty dots represent the mean square dispersion a g 
of the normalized number of groups for TV = 100. 

with our prediction, there is practically no dependence of 
the transition width on the population size TV. 

In contrast with the situation for small a, the growth 
of g in the zone of the transition, a > 2, is due to the dis- 
aggregation of the fully synchronized state into clusters of 
mutually synchronized oscillators. Clusters may have var- 
ious sizes and, consequently, the number of clusters may 
vary between realizations. Also, as we discuss later, os- 
cillators which are mutually synchronized in some of their 
phases need not be synchronized in all of them, so that the 
structure of clustering in a given population is not neces- 
sarily the same for each phase. We call this mixed clus- 
tering structure cross synchronization. The phenomenon 
of cross synchronization reveals a very rich form of self- 
organization, also by comparison with the multi-cultural 
configurations of Axelrod's model. 

Empty dots in Fig. [3] represent the mean square dis- 
persion, over phases and realizations, of the normalized 



number of groups, for TV = 100 and F — 50. We recognize 
the typical fluctuation peak around the critical point, in- 
dicating -in our case- that the number of groups is less 
well-defined in the transition range where synchronization 
breaks down. A more detailed picture of the statistics of 
groups at the transition is provided by FigHl where we 
show histograms of the individual values of the normal- 
ized number of groups for several values of a. As a grows, 
the histogram shifts from a sharp peak at g » to a sharp 
peak at g « 1 , passing by broader distributions of variable 
width as the transition proceeds. 

The phenomenon of cross synchronization addressed to 
above -i.e., the formation of different clustering structures 
for different phases- can be quantified using the normal- 
ized Hamming distance h. In fact, if the distribution of 
oscillators into clusters is the same in all phases, the ma- 
trices m(- (f = l,...,F) are all identical, and the normal- 
ized Hamming distance vanishes. On the other hand, dif- 
ferent clustering structures for different phases generally 
give a positive value of h. For a quantitative evaluation of 
cross synchronization, the value of h must be compared 
with the normalized Hamming distance h r for a random 
distribution of phases. Numerically, h r can be computed 
by randomly shuffling the phases of all oscillators. Fig- 
ure [5] shows, as full dots, the ratio h/h r as a function of 
a for TV = 100 and F = 50 on a random network. Be- 
low the transition, as expected, h/h r — 0. In fact, in the 
fully synchronized state only one group -containing all 
the population- is formed in all phases. For a 2, on the 
other hand, h/h r w 1, showing that well above the tran- 
sition the unsynchronized state is statistically equivalent 
to a random distribution of phases. In the intermediate 
zone, in contrast, h/h r exhibits a sharp maximum attain- 
ing values close to 20. It is therefore in this region that 
cross synchronization is most conspicuous. 
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Fig. 5. Normalized Hamming distance ratio, h/h r , as a func- 
tion of q, for N = 100 and F = 50 on a random network 
(full dots), and on the one-dimensional array, either excluding 
(empty dots) or including (dotted line) twisted states. 

Cross synchronization occurs also on the one-dimen- 
sional array. In this case, however, its detection through 
the Hamming distance is more tricky due to the pres- 
ence of twisted states below the transition. In realizations 
where some of the phases are synchronized while other are 
distributed in twisted states the clustering structures are 
maximally different and, consequently, h attains very large 
values. The dotted line in Fig. [5] represents the ratio h/h r 
as a function of a for realizations on the one-dimensional 
array. To compare with random networks, thus, we exclude 
from the calculation of h the realizations where twisted 
states are found. Empty dots in Fig. [5] stand for the result 
of this procedure. The phenomenon of cross synchroniza- 
tion is qualitatively the same as, but stronger than, in 
random networks. 



4 Conclusion 

Our study of synchronization in a population of coupled 
multi-phase oscillators has been motivated by the affinity 
between the emergence of cultural domains in Axelrod's 
model [1, 2,3] and the occurrence of coherent behaviour in 
distributed interacting systems OHH]. The state of each 
oscillator is characterized by a set of phases, which repre- 
sent the cultural features of Axelrod's model. We focused 
our attention on the long-time asymptotic organization of 
the oscillators into synchronized or unsynchronized config- 
urations, under the effect of pair interactions which are at- 
tractive if the individual states are similar, but are absent 
if the difference between individual states grows beyond a 
given threshold. 

Analytic calculations suggest that, as a function of the 
interaction threshold, the population undergoes a synchro- 
nization transition. If the threshold is too small, individual 
states are generally too different from each other, and the 
population fails to collectively synchronize. Conversely, a 



large threshold enables the attractive interactions and pro- 
motes synchronization. The same calculations show that 
the transition becomes sharper if the number of phases per 
oscillator grows. On the other hand, the transition width 
is independent on the population size. 

These analytical predictions were confirmed by numer- 
ical realizations of our system, both on ordered arrays and 
on random networks. Numerical results also revealed that, 
regarding the organizational structure of the population, 
the most interesting regime occurs precisely around the 
transition. We first found that, in this region, the synchro- 
nized state undergoes disaggregation into clusters of mu- 
tually synchronized oscillators, as observed in other sys- 
tems of interacting oscillators [9] . The number of clusters 
and, consequently, their size, varies with the interaction 
threshold. 

Remarkably, moreover, the segregation of oscillators 
into clusters is not necessarily associated with the syn- 
chronization of all the individual phases. In particular, two 
oscillators may belong to the same synchronized cluster 
with respect to some of their phases -which, in the asymp- 
totic state, adopt identical values- but to different clus- 
ters with respect to the other phases. This phenomenon, 
which we have called cross synchronization, discloses a 
highly complex form of collective organization, consisting 
of coherent clusters mixed with respect to the internal 
variables which specify the individual states. It is inter- 
esting to note that Axelrod's model does not exhibit a 
regime corresponding to cross synchronization. This kind 
of regime would however make sense as a possible distri- 
bution of culture diversity over a real population, in the 
form of partially overlapping cultural domains. We have 
here presented a preliminary quantitative analysis of cross 
synchronization -aiming, mainly, at giving a means to de- 
tect the phenomenon- but, clearly, much further work is 
necessary to fully characterize its statistical features. 
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